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Abstract 

We present Wind spacecraft observations of ion distributions showing field-aligned beams 
(FABs) and large-amplitude magnetic fluctuations composed of a series of shocklets and short 
large-amplitude magnetic structures (SLAMS). We show that the SLAMS are acting like a local 
quasi-perpendicular shock reflecting ions to produce the FABs. Previous FAB observations re- 
ported the source as the quasi-perpendicular bow shock. The SLAMS exhibit a foot-like magnetic 
enhancement with a leading magnetosonic whistler train, consistent with previous observations. 
The FABs are found to have T b ~ 80-850 eV, V,,/V„. ~ 1-2, T^/Tp ~ 1-10, and n,,/n, ~ 0.2- 
14%. Strong ion and electron heating are observed within the series of shocklets and SLAMS 
increasing by factors >5 and >3, respectively. Both the core and halo electron components show 
strong perpendicular heating inside the feature. 



1 Introduction 



Collisionless shock waves are a ubiquitous phenomena in space plasmas and are thought to mediate 



and possibly produce some of the highest energy cosmic rays in the universe \Blandford and Eich- 
|fer||1987) . They are also efficient accelerators of much lower energy particles [e.g. Kis et aZ.||2007[ 
producing a region upstream of the shock transition called a foreshock. The terrestrial foreshock 
has been examined in detail [e.g. \Eastwood et aL 2005 and references therein]. In the process of 
studying this region, five ion populations have been identified: (1) field-aligned beams (FABs), (2) 



intermediate ions, (3) diffuse ions, (4) gyrating ions, and (5) gyrophase-bunched ions [e.g. Kis et al. 
2007 1. All the ion species have been observed to be spatially well separated ^Meziane et al. 201 1) . 



Early observations \Greenstadi\ \1976l suggested that FABs had their origin on field lines con 
nected to a quasi-perpendicular (shock normal angle, or S „, >45°) portion of the bow shock. FABs 
are primarily composed of protons streaming along the ambient magnetic field away from the bow 
shock with temperatures (T,,) ^80-600 eV, densities (n b ) 1-10% of the ambient solar wind density 
(n D ), and beam speeds up to 5 times the solar wind speed (V„) \Bonifazi and Moreno 1981a|b 



Paschmann et al 



1981 1. FABs often have strong temperature anisotropics with T ±b /Tn b > 4-9 
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\Paschmann et al. 1 98 1 1 . They are typically observed in the absence of or near large magnetic 



fluctuations, but not simultaneously with these waves \Kis et al. 2007 Meziane et al. 201 1 

The terrestrial foreshock on the quasi-parallel side has a broad spectrum of large amplitude 
waves. This spectrum consists of transverse Alfvenic waves, right-hand polarized (in plasma frame) 
ultra-low frequency (ULF) waves near the ion gyrofrequency (f ci ), compressional magnetosonic 
waves \ Hoppe and RusseU\\\9%2>\ , magnetosonic-whistler mode waves \Hopp e~et a/.[|1981^ , and an 
ensemble of higher frequency (f > 5-10 Hz) waves up to the electron plasma frequency (f pe ) [e.g. 
\Briand\^2009 and references therein]. We will focus on two specific types of waves in the foreshock, 
shocklets [Hoppe et al. 198 1 1 and short large-amplitude magnetic structures (SLAMS) {Schwartz 



et al. 



1992). Both types are thought to grow out of the ULF wave field, due to an interaction with 

2003 1 . Both can radiate a dispersive higher 



gradients in the diffuse ion densities \Scholer et al. 



frequency electromagnetic whistler precursor due to wave steepening and are always observed si- 

and references therein]. 



2009 



multaneously with diffuse ion distributions [e.g. \ Wilson III et al. 

Though SLAMS and shocklets have different scale sizes ^Lucek et al. 2002 1, they are more eas- 
ily differentiated by their relative amplitudes. Shocklets exhibit weak magnetic compression with 
8B/B„ < 2, while SLAMS are much stronger with 5B/B„ > 2 and often exceeding a factor of 4. 

Previous studies have found FABs near SLAMS, but instrumental limitations prevented an exam- 
ination of the evolution of the ion distributions across the SLAMS [ Schw aTtz etal.\ [l 992 ; Wil kinson] 
\et al.\ [T993) . More importantly, these studies showed indirect evidence that the SLAMS, not the 
bow shock, were producing the FABs locally. 

In this paper we report the first high time resolution observations of the evolution of FABs 
through large amplitude magnetic fluctuations, identified as shocklets and SLAMS, in the terres- 
trial foreshock. The ion beams are more intense on the upstream (sunward) side of the SLAMS, 
suggesting a local source. The paper is organized as follows: Section|2]discusses the data sets and 
analysis techniques, Section [3~T| provides an overview of relevant parameters, Section 3.2 discusses 
the particle distribution observations, and Section |4]presents our conclusions. 



2 Data Sets and Analysis 



The magnetic field was obtained from the Wind dual, triaxial fluxgate magnetometers $Lepping\ 
et al. 1995 1 sampled at ~11 samples/s. Full 4n steradian low energy (<30 keV) ion and electron 
distributions were obtained from the Wind/3DP EES A and PES A particle detectors | Lin et aZ.||1995| 



For more details about the analysis of data from the 3DP instrument, see Wilson III et al. | 2010| . 
The solar wind velocity (V. t „), average ion thermal speed (V77), and average ion temperature (T,) 
were determined with the 3DP PES A Low and SWE Faraday Cups (FCs) IQgilvie et a~L\ [1955). 
Absolute electron densities were determined from the plasma line in the WAVES thermal noise 



receiver (TNR) instrument Bougeret et al. | 1995[ and used as an estimate of the solar wind ion 
density (n,). To analyze the beams, we fit the FABs to bi-Maxwellians to determine estimates of the 
density (n fc ), the temperature (T b ), and beam speed (V;,). We will use these values to compare to 
previous observations and provide better constraints for theory. 

The wave vector, k, and the polarization with respect to the quasi-static magnetic field were 



determined using Minimum Variance Analysis (MVA) \ Khrabrov and Sonnerup 1998 1. The details 
of this technique are discussed in \Wilson III et al. | 2009) . We calculated the angles between the 
wave vector and the local magnetic field (9 kB ) and the solar wind velocity (9 kv ). 

We used two methods to determine the bow shock normal vector. The first involved the use of 



the Rankine-Hugoniot conservation relations \Koval and Sza bo 1 2008 1 with the parameters observed 
at the last crossing of the bow shock. From this we derive a single shock normal vector. The second 
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method involved projecting the local magnetic field vector onto the surface of a model bow shock 
\Slavin and Holzer\ \[9%\ |. Once we determined the shock normal vector, we were able to calculate 
the shock normal angle, 9 B „. 



3 Observations 



3.1 Foreshock Observation Overview 

Figure [T] plots an overview of the magnetic field measurements observed for the two events of in- 
terest, which show structures identified as groups of shocklets and SLAMS. We will focus on the 
2002-08-10 event herein. First we needed to determine whether the spacecraft was magnetically 
connected to the terrestrial bowshock and then eliminate the possibility that these structures were 
simply due to an expansion of or close encounter with the terrestrial bow shock. Wind was located 
at a GSE position of ~<+13.6, -12.4, +0.02> R £ for the group of SLAMS in Figure [T]C-E and 
at ^<+13.8, -3.0, +0.4> R £ for Figure[TJl-J. These positions correspond to distances of ^2.4 R E 
and ~1.4R £ from the last bow shock crossings, respectively. 

To exclude the possibility of a bow shock expansion, we examined data from the ACE, GOES 
8 and 10, and Interball spacecraft (courtesy of CDAWeb). We found no transient features in any of 
these data sets suggesting a sudden reduction in solar wind pressure occurred that could allow the 
bow shock to expand >1 R £ out to the Wind position at the time of interest in either event. There- 
fore, we conclude that these structures are features of the terrestrial foreshock and not an expansion 
of the bow shock. 

Next, we needed to determine the local geometry of the bow shock. The blue line in Figures 
[TJi and[TJj represents the 9 B „ calculated from our first method and the magenta line was calcu- 
lated using the second method, both discussed in Section [2] Note that the second method used a 
smoothed average magnetic field, instead of the HTR MFI data, and gaps indicate regions not mag- 
netically connected to the model shock surface. The Rankine-Hugoniot solutions gave us shock 
normal vectors of -<+0.750, -0.621, -0.134> for the 2000-04-10 event and -<+0.987, -0.158, 
+0.002 > for the 2002-08-10 event. We estimated the upstream foreshock average magnetic field, 
B,„ using the Wind/MFI observations between 14:20-15:22 UT for the 2000-04-10 event and be- 
tween 11:21-12:49 UT for the 2002-08-10 event. The average GSE vectors were found to be B„ ~ 
(+3.60,-2.99,-1.47) nT and ~ (-4.32, +0.65,-0.41) nT, respectively. These estimates with the 
above shock normal vectors give 9 B „ ~ 14° and ~ 6°, respectively. Therefore, the Wind spacecraft 
was primarily traveling through the quasi-parallel region of the terrestrial foreshock for the time of 
interest. 

The foreshock structures, marked by vertical red lines in Figures [T]\-B and [TJ^-G, were com- 
posed of a series of compressive magnetosonic waves (5B in phase with density, 8n, fluctuations) 
identified as shocklets and SLAMS shown in Figures [T]C-D and [TJl-I, respectively. The SLAMS 
were observed to have: (1) mixtures of right- and left-hand polarizations (spacecraft frame); (2) 
very oblique (9 kB > 55° and 9 kv > 40°) propagation; and (3) 5B/B„ > 2-6, consistent with previous 
observations ^Schwartz et~al\ [T9921 \Wilkinson et al\ [T993) . The 2002-08-10 event (Figure [Tp-I) 
shows an isolated SLAMS near 12:52:20 UT The 2000-04-10 event (Figure [T]C-D) did not show a 
similar structure. The importance of this difference will be discussed in the next section. 

Both groups of SLAMS have higher frequency fluctuations on their leading/upstream (i.e. the 
right-hand side of Figures [T]C-D and[Tp-I) edges, consistent with whistler mode waves. The charac- 
teristics of the waves immediately upstream of the the steepened edges are consistent with previous 



observations of whistler precursors [e.g. Wilson III et al. 2009 1. The whistler amplitudes and beam 
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intensity decrease away from the leading edge of the group of SLAMS and eventually the whistlers 
disappear when the FABs disappear. However, we cannot definitively show that the two phenomena 
are causally related because whistler modes have been observed in the absence of FABs. While 



simulations have found that reflected ions can provide free energy for whistler precursors [e.g. Sc 



holer et al. 2003), supported by recent observations [Wilson III et al. 2012 1, their primary source 



is thought to be dispersive radiation [e.g. Sundkvist et al. | [2012) . The source of the whistler mode 
waves are beyond the scope of this paper and we will not discuss them further. 



3.2 Particle Distributions 

We examined the effects caused by the series of shocklets and SLAMS on the ion and electron dis- 
tribution functions for the two foreshock passes. The group of SLAMS created a rarefaction region 
behind the structures with a strong deflection of the solar wind core, analogous to the wake created 
by an obstacle in a fluid flow. The SLAMS caused strong anisotropic heating in the low energy 
(<1.1 keV) electrons and ions (<10 keV). 

Figure[2]shows the HTR MFI data and select PES A High distribution functions (in the solar wind 
frame) for the time range corresponding to Figure [TB-I. The solar wind core is clearly identified in 
the center of panels S-AE and separated from the FABs seen near ^500-900 km/s in panels S-AB, 
moving anti-parallel to B D , which corresponds to the sunward direction for this event. Note that 
panel B contains a FAB as well, but at lower speeds, which can be seen by the difference between 
the foreshock observations in panel A. The FABs on the downstream(earthward) side of the group of 
SLAMS is weaker for both events. Hot diffuse ions (i.e. nonthermal tail observed above ^800 km/s) 
are observed continuously between 12:50:13 UT and 12:51:45 UT, simultaneous with the SLAMS. 
The ion core experienced strong enough heating in this region to be observed by PESA High in pan- 
els D-R in Figure [2] which was supported by the PESA Low and SWE distributions. These effects 
are consistent with previous observations | Schwartz et al. 1992 Wilkinson et al. 1993 1. 



The ion beams shown in Figure [2] are intense FABs observed between 12:51:51-12:52:48 UT 
(At FAB ~ 57s), which is a considerably longer duration than for the 2000-04-10 event (At^ ~ 10s) 
but comparable to previous observations [e.g. Meziane et al. 2004) . The FABs between 12:51:51— 
12:52:48 UT have T b ~ 175-850 eV, T ± JT U > 2.3-9.7, V b /V„ ~ 1.1-2.4, and n,,/n, - 0.3-14%, 
consistent with previous observations [e.g. Bale et al. 2005 and references therein]. Correspond- 
ingly, the FABs observed at the 2000-04-10 event (not shown) had T A ~ 80-200 eV, V,/V. nl , ~ 
1.3-2.0, and n,/n, ~ 0.2-1.6%. Though we did examine theoretical growth rates associated with 



ion/ion beam instabilities [e.g. Akimoto et al. 1993 1, we do not believe the FABs are the source of 



the SLAMS, rather we will show that the SLAMS are the source of the FABs. We believe the source 
of free energy for SLAMS growth are the diffuse ions observed simultaneously with each group 



of SLAMS, consistent with theory [e.g. Scholer et al. 2003 1 and previous observations {Schwartz 
|eFan|l992) . 

Note that FABs were observed both downstream (earthward) and upstream (sunward) of the 
group of SLAMS for both events. The FAB intensity was found to be greater on the upstream(sunward) 
side of the group of SLAMS but weaker on the upstream side of the isolated SLAMS (^12:52:20 
UT) shown in Figure[Tp-I. This suggests that the group of SLAMS may act like a quasi-perpendicular 
shock, shown in Figures [Tp andflj, consistent with previous observations [M ann et q/.||1994) . The 
decrease in FAB intensity upstream of the isolated SLAMS suggests that groups are necessary to 
create an obstacle strong enough to reflect a significant population of ions, consistent with theory 
[Schwartz and Burgess 1991]. If the FABs were produced at the bow shock, then it is very surpris- 
ing that the beams are so coherent after traversing the turbulent groups of SLAMS. If the magnetic 
field cannot connect the spacecraft with the bow shock without going through the SLAMS, then the 
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SLAMS must be the source of the FABs. 

Figure [3] is an illustrative cartoon that we will use to argue that the group of SLAMS in each 
event is the source of the FABs. If we assume the group of SLAMS, referred to as the obstacle 
for brevity, are being convected with the solar wind at roughly V w , then the spacecraft will be ef- 
fectively stationary with respect to the obstacle. Therefore, the path of the spacecraft through the 
obstacle is anti-parallel to V.„„ and of length L, = V„, At sc . The spacecraft will be shielded from the 
terrestrial bow shock along magnetic field lines for a distance L shadow , Since we know that V.„, is at 
an angle to the background magnetic field, B„, then L shadow >(=) L s for angles <(=) 90°. It takes 
AU — 136 s(108 s) to traverse the obstacle in the 2000-04-10(2002-08-10) event, much greater than 
the corresponding At FAB discussed above. The average complementary angle between V„ and B„ is 
< 45° for both events, therefore L a > L s . This means that the spacecraft is in "magnetic shadow" 
of the obstacle for at least At sc after exiting on the upstream side. Therefore, the spacecraft cannot 
connect to the bow shock without going through the obstacle. The amount of turbulence and rota- 
tion in the magnetic field observed through the obstacle suggests that a bow shock source would be 
highly unlikely. In conclusion, we argue that the group of SLAMS must be the source of the FABs 
observed on their upstream side, not the bow shock, consistent with predictions of the quasi-parallel 
bow shock ^Schwartz and Burgess\\199l\ . 



4 Discussion and Conclusions 



This study presents observations of the evolution of field-aligned ion beams (FABs) through SLAMS 
in the terrestrial foreshock. The FAB intensities were higher on the upstream(sunward)-side of the 
group of SLAMS than the downstream. The group of SLAMS created an effective wall between 
the spacecraft and the bow shock, which suggests that these structures are locally producing the ob- 
served beams. However, as we discussed the isolated SLAMS near 12:52:20 UT in the 2002-08-10 
event may be too thin to locally produce FABs. In addition, the field is observed to rotate to a quasi- 
perpendicular orientation within the SLAMS, which would support the previous statement and the 



predictions of \Schwartz and Burgess] fl991) . 

The FABs propagate (in the plasma frame) away from the bow shock, consistent with previous 
observations [e.g. Kis et al. 2007 1 . This is also the direction toward the upstream(sunward)-side 
of the SLAMS. They had T b ~ 80-850 eV, V fc /V„ - 1.1-2.4, and n fc /n, - 0.3-14%, consistent with 
previous observations [e.g. Bale et al. 2005 and references therein]. While ion beams have been 



previously observed near SLAMS and suggested to be locally produced \Schwartz et al. 1992 



Wilkinson et al. 1993[ , no previous reports have shown the evolution of FABs through SLAMS and 



showed the SLAMS to be the source. The group of SLAMS was shown to create a barrier between 
the spacecraft and the bow shock along magnetic field lines, which supports our argument of local 



beam production and confirms predictions of Schwartz and Burgess \ 1991) . 

The source of the SLAMS and their associated whistler precursors is beyond the scope of this 
manuscript, however simulations suggest that diffuse ions are responsible for SLAMS and the com- 
bination of dispersion and reflected ion-driven instabilities explain the precursors \Scholer et al.\ 
2003; Sundkvist et al., 2012 1. The simultaneous observation of diffuse ions with the SLAMS and 
FABs with the precursors does support these results, but precursors have been observed in the ab- 
sence of reflected ions. 

In conclusion, we show the first direct evidence that groups of SLAMS can act like a quasi- 
perpendicular shock producing reflected field-aligned ion beams. 
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Figure 1 : The top half of the figure shows three second resolution of the magnitude and the normal 
incidence frame (NIF) Sundkx nst et al.\ [e.g. 2012 1 components of magnetic field data from the Wind 
spacecraft on 2000-04-10 (A-B) and 2002-08-10 (F-G) each with three vertical lines that indicate the 
magnetopause crossing (green), the last bow shock crossing (blue), the foreshock boundary (purple), 
and the red lines show the time periods for panels C-E and H-J. The tick mark labels at the bottom 
of these two panels are: UT time, and the Wind spacecraft radial distance (Re), GSE longitude 
(degrees), and GSE latitude (degrees). Every panel has the same format, but the bottom two panels 
show the HTR MFI data and the shock normal angle for model (magenta) and Rankine-Hugoniot 
(blue) solutions. 
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Figure 2: Selected PESA High Burst distributions shown for the time range shown in Figure [Tp- 
J. The ion distribution plots are contours of constant phase space density (uniformly scaled from 
lxl(T 13 to lxl(T 9 s 3 cm 3 km 3 , where red is high) with arrows showing projections of V,,, 
(black), sun direction (blue), and Earth direction (magenta). Panels S-AE have circles of constant 
energy at 500, 700, 900, and 1100 km/s. The coordinate system used for the contours is shown in 
the lower right-hand corner of the plot, where each plot ranges from ±1500 km/s on each axis. Ion 
beams are clearly identified in panels B and S-AB. 
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Figure 3: A schematic cartoon used to illustrate how the series of SLAMS in Figures [T|C-E and 
[Tp-J can block the spacecraft (SC) from "seeing" the bow shock along magnetic field lines. In this 
example, the sun is to the right and Earth to the left. 



11 



